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Abstract 

The Standard-Model Higgs boson with mass M# 3> 2Mz decays almost exclu- 
sively to pairs of W and Z bosons. We calculate the dominant two-loop corrections, 
of O(GpMfj), to the partial widths of these decays. In the on-mass-shell renor- 
malization scheme, the correction factor is found to be 1 + 14.6%(M#/TeV) 2 + 
16.9%(M#/TeV) 4 , where the second term is the one-loop correction. We give full 
analytic results for all divergent two-loop Feynman diagrams. A subset of finite 
two- loop vertex diagrams is computed to high precision using numerical techniques. 
We find agreement with a previous numerical analysis. The above correction factor 
is also in line with a recent lattice calculation. 
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1 Introduction 



One of the longstanding questions of elementary particle physics is whether nature makes 
use of the Higgs mechanism of spontaneous symmetry breaking to endow the particles 
with their masses. The Higgs boson, H, is the missing link sought to verify this theoretical 
concept in the Standard Model (SM). Its mass, Mh, is essentially a free parameter of the 
SM. The failure of experiments at the CERN Large Electron-Positron Collider (LEP 1) 
and the SLAC Linear Collider (SLC) to observe the decay Z — > ffH has ruled out the 
mass range Mh < 65.2 GeV at the 95% confidence level jl]. The discovery potential of 
LEP 2, with centre-of-mass (cm.) energy yfs = 192 GeV and luminosity L = 150 pb _1 per 
experiment, will reach up to M H « 95 GeV JJ. We note in passing that M H lower bounds 
may also be derived theoretically by demanding that the Higgs vacuum be (meta) stable 

!• 

Several theoretical arguments bound Mh from above. The requirement that partial- 
wave unitarity in intermediate-boson scattering at high energies be satisfied at tree level 
establishes an upper bound on M H at about (8ny/2/3G F ) 1/2 « 1 TeV ||. This unitarity 
bound is significantly lowered by including one- || and two-loop || corrections. However, 
the improved bound depends on the energy scale, A, up to which the SM is assumed to 
remain valid. The same is true for the triviality bound, which is derived perturbatively 
by requiring that the running Higgs self-coupling, A(/i), stays finite for renormalization 
scales fi < A. Lattice computations, which take all orders of the perturbative expansion 
into account, yield a triviality bound of M H ~ 700 GeV M. The precise value of this 
bound depends on the specific form of the lattice action and the regularization method 
used. 

An alternative way of constraining Mh from above is to require that the Higgs sector 
be weakly interacting, so that perturbation theory is meaningful |§ . The resulting bounds 
depend somewhat on the considered process and the precise definition of perturbation- 
theory breakdown, but they are independent of assumptions concerning the scale A beyond 
which new physics operates. At one and two loops, the leading high-M^ corrections to 
physical observables related to Higgs-boson production or decay are of 0(G F M H ) and 
0(G 2 F M H ), respectively. The one-loop corrections to the partial widths of the Higgs- 
boson decays to pairs of fermions [1(J and intermediate bosons |ll| are relatively modest, 



below 15% at Mh = 1 TeV if the on-mass-shell renormalization scheme is employed. The 
situation is similar for Higgs-boson production |12|]. For Mh increasing, the 0(G F M H ) 



corrections will eventually exceed the 0(G F M H ) ones in size. In the case of the fermionic 
Higgs-boson decays, this happens at Mh ~ 1.1 TeV in the on-shell scheme [^,|T^]. Obvi- 
ously, the perturbative expansion in G F M H ceases to usefully converge for larger values 
of Mh- The study of high-energy weak-boson scattering at two loops (with neglect of 
trilinear couplings) leads to a perturbative upper bound on M H which depends on the 
considered cm. energy. For y/s in the 1-TeV range, the Mh upper bound comes out as 
low as 450 GeV [U| . However, this result may be relaxed using resummation techniques 

I- 

Since a Higgs boson with Mh 3> 2Mz decays almost exclusively to W + W and ZZ 
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pairs, it is desirable, both from the theoretical and phenomenological points of view, to 
find the two-loop 0{G 2 F Mfj) corrections to the partial widths of these decays as well. This 
is the purpose of the present paper. In contrast to the fermionic case JT^JT^] , this involves 
the calculation of massive three-point functions at two loops. Due to the presence of a 
large external momentum, heavy-mass expansion techniques are not applicable here. As 
in a previous analysis ||13|| , we take advantage of the Goldstone-boson equivalence theorem 
17| and neglect the gauge and Yukawa couplings as well as the intermediate-boson masses. 



In this approximation, the radiative corrections to the H — > W + W~ and H — > ZZ decay 
widths coincide. 

In the evaluation of the relevant two-loop Feynman diagrams, both infrared (IR) and 
ultraviolet (UV) divergences occur. We calculate all divergent diagrams analytically, us- 
ing dimensional regularization. This guarantees that the cancellation of the IR and UV 
divergences is manifest. The remaining diagrams are individually devoid of singularities. 
They are computed numerically in four space-time dimensions after formulating the inte- 
grals in parallel and orthogonal space. This has the advantage that potential numerical 
instabilities in connection with the conventional Feynman parameterizations are avoided. 
We list our results for the individual diagrams, so that they are available to other authors 
for comparisons and further calculations. 

We anticipate that we find good agreement with a recent, pioneering work by Ghin- 



culov ||18|| , which fully resorts to numerical methods. In Ref. |[L8|| , the IR divergences are 
regularized by introducing a small Goldstone-boson mass. The final result is obtained 

is 



by numerically taking the massless limit. Another conceptual difficulty in Ref. [|18[ 
related to the treatment of (endpoint) singularities in the Feynman-parameter integrals. 
These are controlled by choosing suitable integration paths in the complex plane with 
the aid of spline functions. The values of the individual diagrams are not specified in 
Ref. |lBfl . We believe that, in view of the enormous complexity of the problem at hand, 



it is indispensable to have at least two independent analyses that agree with each other. 

This paper is organized as follows. In Sec. 2, we catalogue the contributing Feynman- 
graph topologies and explain our renormalization framework. In Sec. 3, we describe our 
analytical and numerical techniques. In Sec. 4, we present our results for the individual 
two-loop vertex diagrams and assemble the renormalized H — > VV (V = W, Z) transition 
amplitude. In Sec. 5, we explore the phenomenological implications of our result and draw 
the conclusions. The relevant one-loop three-point amplitudes to 0(e), where D = 4 — 2e 
is the dimensionality of space time, are listed in the Appendix. 



2 Framework 

In this section, we set the stage for our calculation of the H — > VV (V = W, Z) de- 
cay widths to 0{G 2 F Mf I ). As mentioned in the Introduction, this task may be greatly 
simplified in the limit of interest, Mh 3> 2Mz, through the use of the Goldstone-boson 
equivalence theorem jnj. This theorem states that the leading high-M# electroweak con- 
tribution to a Feynman diagram may be calculated by replacing the intermediate bosons 
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W ± , Z with the would-be Goldstone bosons w , z of the symmetry-breaking sector of the 
theory. In this limit, the gauge and Yukawa couplings may be neglected against the Higgs 
self-coupling. By the same token, the Goldstone bosons may be taken to be massless, and 
the fermion loops may be omitted. 



Adopting the conventions of Ref. 



we may write the relevant Lagrangian for the 
symmetry-breaking sector of the SM in terms of bare quantities, marked by the subscript 
0, as 



C 



SBS 



~<9 M w ■ 0"w o + ^HoVHo - \m 2 



u),0 W 



2 - -M 2 Hfi H 2 



Ao 
4 



(w 2 + Hi) 2 - \ v (w 2 + Hi) H 



where A is the Higgs quartic coupling, v is the Higgs vacuum expectation value, H is 
the real scalar Higgs field, and the real scalar SO(3) triplet, w = (wi, w^, w 3 ), is related 
to the Goldstone bosons, w ± and z, by w ± = (wi =F iw 2 )/V^ and z = w 3 , respectively. 
The tadpole counterterm, which cancels all tadpole contributions of £q BS order by order, 
has been omitted in writing Eq. ([I]). Consequently, all diagrams which include tadpole 
contributions need to be dropped in calculations |19| . 

We work in the on-mass-shell renormalization scheme as formulated in Refs. |13,19[ . 
The residual SO (3) symmetry present in Eq. ([I]) is preserved in the renormalization pro- 
cedure, so that the w and z bosons share one mass counterterm and one wave-function 
renormalization constant. Specifically, the renormalization conditions are as follows: (i) 
The mass counterterms, SM^ and SMfj, are determined so that the renormalized masses, 
M 2 = M 2 — 5M 2 and Mjj = Mfj — 5M]±, coincide with the physical (pole) masses, i.e., 
the pole positions of the radiatively corrected propagators. In the case of the Goldstone 
bosons, one has M 2 = 0. (ii) The wave-function renormalization constants, Z w and Zh, 
which enter the relations = Z}J 2 w ± and H = Z\j 2 H between the bare and renor- 
malized fields are adjusted so that the on-shell propagators have unit residues. 



(iii) The 

physical Higgs quartic coupling and vacuum expectation value are fixed by requiring that 
A = GpMfj/y/2 and v = 2~ 1//4 G F 1//2 , respectively, where Gf is Fermi's constant. This 
leads to the following relations jl3l|l9|: 



5M 2 = 


-n°(o), 






5M 2 H = 


-ReH° H (il4), 






1 

Z w 


i _ A n ° (p 2 ) 
dp 2 wyp ' 




o' 


1 

z^ = 


1 - ±ReU W) 


p 2 =Mjj 


A = 


A / 6M 2 H -SM 2 \ 
Z w \ + Mjj ) 


v = 









(2) 



where n° (p 2 ) and U° H (p 2 ) are the self-energy functions of the bare Goldstone and Higgs 
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fields, respectively, calculated from Lagrangian (|1|). Two- loop expressions for SM^, 5Mfj, 
Z w , Zh, and Aq in terms of Mh and A may be found in Refs. [|T9| , pO| . 



The Feynman diagrams that contribute to the H — » VV transition amplitude through 
0{G%Mjj) are depicted in Fig. [l|. Dashed and solid lines represent Goldstone and Higgs 
bosons, respectively. Adjacent propagators with identical four-momenta, which occur in 
connection with mass-counterterm insertions, are separated by solid circles. The calli- 
graphic labels stand for the values of the diagrams and are used in our formulas; the 
nomenclature extents the one introduced in Ref. [19]. The underlying conventions are 
explained in Sec. 3 and the Appendix. 

There are 6 one-loop diagrams (<2, £>j, and Cj), 18 reducible two-loop diagrams {BiBj, 
BiCj, QTi, and QT>i), and 52 irreducible two-loop diagrams (Ai, £«, Li, Ti, Qi, Hi, 
Zj, J7i, and /Cj). In addition, there are 8 one-loop diagrams with one mass-counterterm 
insertion {5Ml% and 5M%V t , where a = w,H). Further reducible two-loop diagrams 
emerge as products of one-loop vertex and wave-function-renormalization diagrams; these 
are not counted here. The one-loop diagrams need to be evaluated through 0(e). The 
corresponding results for the tadpole (Q) and two-point topologies (£>; and %) may be 



found in Ref. [19]; those for the three-point topologies (Cj and T>i) are listed in the 



Appendix. The irreducible two-loop two-point topologies (Ai, Si, and £j) may be taken 



from Ref. fl9|| . Here, we regard two external legs that merge into the same vertex as one 
point; solid circles are not counted as points. For example, we view both % and Ai as 
two-point topologies. 

The main technical difficulty of the present paper is to tackle the irreducible two-loop 
three-point topologies (Ti, Qi, Hi, Xj, Ji, and /Q). Topologies Ti, Qi, and Hi are ultraviolet 
divergent. Topologies H3 and H4 are also plagued by infrared singularities. Topologies 
Ji, Ji, and /Cj are UV and IR finite. After describing our computational techniques in 
Sec. 3, we itemize the results for topologies T\ through K 2 in Sec. 4. 

Calculating the combinatorial factors with the help of the software package DIA- 
GRAMMAR |nj, we find the following master formula for the two- loop H — > VV tran- 
sition amplitude in terms of A , v , and M H : 

T(H -> VV) = -2Xov Z w z]l 2 jl + i\ (10B + 8B 1 + QB 2 ) + i(A « ) 2 (8C 1 + 24C 2 ) 

- Xl(50B B + 3QB B 2 + + 18B 2 B 2 + 20QT + 8QTx + 24QT 2 + 36QT 3 

+ 4QA 10 + 4(U im + 16^20 + 12^ 2m + 24^3 + 36A 4 
+ 60^1 + 4JF 2 + 24J-3 + 40^4 + 16JS + 24^ 6 ) 

- Ao(A vo) 2 (40i3 C 1 + 2AB C 2 + 32B 1 C 1 + 96SiC 2 + 2AB 2 d + 12B 2 C 2 

+ 1QQV 1 + 2AQV 2a + 2AQV 2b + IUQV 3 

+ 80£i + 32£ 2 + 48£ 2 * + 72£ 3 + 144£ 4 + 216£ 5 

+ 40£i + 32£ 20 + 144£ 2m + 96£ 3 + 216£ 5 

+ 288^1 + 288£ 2 + 64£ 3 + 96£ 4 + 48£ 5 + 80Q 6 

+ 48Ji + 48J 2 + 32J 3 + 288J 4 + I6OJ5 + 96X 6 + 32X 7 ) 

- (A w ) 4 (864H 1 + 288^2 + 64ft 3 + 96H 4 + IUH 5 + A8H 6 + 288/Ci + 32/C 2 
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+ 



+ 
+ 



864 Ji + 96 J 2 + 96 J 3 + 32 J" 4 + 576 J 5 + 192 J % + 192 J 7 + 64J" 8 ) 
+ 87^) + i(\nVn?(16Vi + 2AVo h )} 5M 2 , 



iA (20T + 8Ti) + i(A t;o) 2 (16Pi + 24X? 26 )" 
:A (8T 2 + 12T 3 ) + z(A ^o) 2 (8^ 2 a + 48X? 3 )] 5M£} . (3) 



Here, the tree-level term has been factored out together with the wave-function renor- 
malization constants of the external legs; in our two-loop approximation, these quantities 
need to be evaluated through 0(A 2 ). The second and third terms within the curly bracket 
arise at one loop; their prefactors of Ao and Vq need only to be expanded through 0(A). 
The residual terms are of genuine two-loop order; here, Ao and Vq may immediately be 
replaced with A and v, and it is sufficient to insert the one-loop expressions for SM 2 
and SMfj. Notice that we have already performed mass renormalization in Eq. (|3|), so 
that Mh, which appears inside the expressions for the various diagrams, is the only mass 
parameter left. 



3 Computational techniques 

As usual, we calculate the Feynman diagrams in D = 4 — 2e space-time dimensions 
and introduce an unphysical 't Hooft mass scale, /i, to keep the Higgs quartic coupling 
dimensionless. After coupling renormalization, Eq. (|3p will be finite and fi independent 
in the physical limit e — > 0. To fix the notation, we write down the generic one- loop 
two-point integral (with two scalar propagators): 

B(p\ ml ml) = /i 2e J ^ q2 _ m 2 + ie) ^ q + v)2 _ m 2 + ie] - (4) 

The one-loop two-point diagrams in Fig. p] are then defined as 

B = B(M 2 H , 0,0), 
B 1 = B(0,0,M 2 H ), 
B 2 = B(M 2 H ,M 2 H ,M 2 H ), 

^ = ^f{M 2 H ,0,0), 

ri = ^B(o,o,0, 

d 



d 



°2 

T 3 = -^B(M 2 H ,M 2 H ,M 2 H ). (5) 



Analytic results for these diagrams to 0(e) may be found in Ref. [19]. Note that, in 



contrast to Ref. []nj , we exclude the factors of i connected with the scalar propagators in 
the definitions of the diagrams. For the one- loop tadpole diagram, we have Q = M\B\. 
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The one- loop three-point diagrams, Ci and Vi, are defined and evaluated to 0(e) in the 
Appendix. 

The two-loop diagrams are defined in analogy to Eqs. (Q) and (p2|) . We first consider 
the IR-finite two-loop diagrams in Fig. |l|. In the limit e — > 0, any member, M, of this 
class may be cast into the generic form 



M 



47ifi 2 \ 2e [A B 



(4tt) 4 (M2)^-4 \M 2 H e\ 



+ C + 0(e) 
e 2 e 



(6) 



where 7 is the Euler-Mascheroni constant, N is the number of propagators in A4, and A, 
B, and C are finite, complex numbers. The special form of the prefactor is to suppress 
the appearance of the familiar term ln(47r) — 7 in B and C . (By contrast, the 7 piece of 
the prefactor was expanded in Ref. Coefficients A, B, and C for diagrams Ai, £i, 

and Li are available in analytic form [[HJ. Diagrams £1 and £2 are IR divergent and will 
be discussed further below. 

The irreducible two-loop three-point diagrams with four propagators, J-'i, exhibit UV 
divergences in the form of quadratic and linear poles in e, but they are IR safe. They may 
be evaluated analytically from appropriate Feynman-parameter representations, proceed- 



ing along the lines of Ref. [19|. In the case of diagram T\, where all internal lines are 



massless, we find agreement with Ref. p2 



Diagrams Gi and Hi are more involved, and it is advantageous to apply dispersion 



relations in a way similar to Ref. P^fl . Let us first concentrate on diagrams Qi. We 
observe that they contain a one-loop two-point subdiagram consisting of a bubble and an 
adjacent propagator. Calling the masses inside this bubble m\ and m,2, the mass of the 
adjacent leg 777,3, and the common four-momentum q, we may rewrite this subdiagram as 

B(q 2 ,m\,ml) _B(ml,m\,ml) i r°° ds y/\(s,ml,m$) 1 + 0(e) 
q 2 — m 2 + ie q 2 — m 2 + ie (An) 2 i(mi+m 2 ) 2 s s — ml — ie q 2 — s + ie 

(7) 

where A(s, m\ , m 2 ,) = [s — (mi + m 2 ) 2 } [s — (mi — m 2 ) 2 ] is the Kallen function. If Eq. (|7|) is 
inserted in the expression for the one-loop seed diagram, the first term turns into a product 
of the type BiCj, which contains all divergences. In the second term, we may interchange 
the dispersion and loop integrations and are left with a finite dispersion integral, which 
may be solved analytically. 

Diagrams Tii may be treated in a similar fashion. They contain a one-loop two-point 
subdiagram consisting of a bubble and two identical adjacent propagators. A useful 
representation of this subdiagram emerges by differentiating Eq. (0) with respect to m\. 
The divergences of diagrams Hi are contained in products of the form B{Dj. In particular, 
this method allows us to separate the IR divergences of diagrams H3 and H4. By the 
way, this method may also be applied to the IR-divergent diagrams £1 and £2- 

Let us now turn to the cancellation of the IR divergences in Eq. @. It is possible 
to extract the IR divergences along with the UV ones as poles in e using dimensional 
regularization. In fact, this avenue was taken in Ref. ||19|| . Here, we adopt an alternative 



strategy. We identify a basic set of IR-divergent one-loop diagrams and show that their 
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prefactors vanish. In this way, we do not need to solve any IR-divergent loop integral. It 
is easy to see that the diagrams which emerge by duplicating any of the Goldstone-boson 
propagators in diagrams Bo, B\, Ci, and C 2 are IR divergent. The resulting diagrams are 
To, Ti, T>i, and V 2 b, respectively. As for the irreducible two-loop diagrams, £\, £ 2 , 
and are IR divergent. Hence there are twelve IR-divergent terms in Eq. (|3|). By means 
of the dispersion-relation method described above, we find that E\ — Bi%, £ 2 — B{T\, 
H3 — B1V1, and 7i 4 — B{D 2b are finite. Exploiting the one-loop identity jl9| 



16M 2 = XQ + ^XvfB^ 



(8) 



it follows on that the twelve IR-divergent terms in Eq. (|3|) may be grouped in four finite 
sets. Specifically, we have 



{\Q-i5M 2 w )T Q + A{\v) 2 £ l 
(XQ~t5Ml)T 1 + 4(Xv) 2 S 2 

(XQ - i5M 2 )V 1 + A(Xv) 2 H 3 = \ 

v 

(XQ-i5M 2 w )V 2b + A(Xv) 2 H A = ^ 



X [-2C(2) + 1 
A [-2 + 0(e)] 



111 



0(e)] 



~C(2)-~ + *f ln2 + 0( e ; 



C(2) v/3 



12 



7T 



2V3 



1 

2 



0(e) 



(9) 



where ( is Riemann's zeta function, with values ((2) = 7r 2 /6 and ((3) ~ 1.202 056 903, 
and Cl 2 is Clausen's integral, with value Cl 2 (7r/3) w 1.014 941606. 

In the remainder of this section, we describe the evaluation of the finite diagrams X\ 
though ¥L 2 . This is achieved with the help of the methods developed in Refs. p4]j25|j26|j27| . 
These methods suggest a separation of the integration space in an orthogonal and a parallel 
space. The parallel space is defined to be the linear span of the exterior momenta involved 
in the process under consideration. For a three-point function, this is a two-dimensional 
space. Each loop momentum is a sum of two vectors, one being its projection onto the 
two-dimensional plane spanned by two non-degenerate exterior momenta, the other one 
being orthogonal to this plane. This splitting does not break Lorentz invariance [2S|. In 
the following, we define orthogonal- and parallel-space variables for our purposes. We 
assume two independent exterior momenta, p\ and p 2 , and two loop momenta, I and k. 
Without loss of generality, we assume p± and p 2 to be timelike. Then, t\ = pi/y'pf and 

&2 = (P2 — £>2-eiei)/ 1 \J (p 2 .e 1 ) 2 — p\ are orthogonal unit vectors, with e\ = —e 2 . = 1 and 
e\.e 2 = 0. We further define 



l = Lei, 
k = k.ei, 
h = l.e 2 , 
ki = k.e 2 , 

l± = l- ke x - he 2 , 
k± = k — koei — k\e 2 . 



(10) 
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We then have 



eji_L = 0, 
ei.k\ = 0, 



72 _ ,2 _ ,2 _ ,2 
i — l O l l '±' 

/c — /lq /c^ A/ 1 , 
Lpi = ZoPi-ei, 
fc.pi = k p 1 .e 1 , 
l.p 2 = l P2-ei ~ kP2-e 2 , 
k.p 2 = k p 2 .e 1 - hp 2 .e 2 , 
l.k = l k — liki — l±.k±, 
(I + kf = (l + hf - (h + h) 2 -ll-k 2 ± - 2l ± .k ± . (11) 

Futhermore, we call the angle between the two-loop momenta in orthogonal space 

(12) 

n\k\ 

In its most straightforward application, this approach delivers a two-fold integral rep- 



resentation for the two-loop two-point function | [25|| . Typically, the remaining integrations 
are with respect to the parallel-space component of each of the two loop momenta. 

For two-loop three-point functions f2l| , and for four-point functions as well PD| , this 
approach opens a route to three-fold integral representations, regardless of the internal 
topology of the function. In Ref. [27(1, this was improved to give a two- fold integral 



representation for the planar three-point case. Remarkably, the so-obtained integral rep- 
resentations are solely defined on finite domains of integration ||27|| . These results are 
sufficient to treat topologies J7i and, with slight modifications, topologies Tj as well. A 
similar result for the three-point non-planar case was obtained recently |[30|| , and was used 
here for the first time to handle topologies /Q. 

We now wish to describe some features of the method. Following separation into 
parallel-space and orthogonal-space variables, the angular integration in the orthogonal 
space is complicated due to the presence of propagators involving both loop momenta, 
which results in a non-trivial integration over the variable z. In fact, the z integration 
delivers a result of the form f(Pi \i ± =k ± =o), where / is a function having a single branch- 
cut and Pi is such a propagator involving both loop momenta. In /, this propagator 
is evaluated at nullified orthogonal loop momenta. Typically, in four- dimensional field 
theory, / turns out to be a logarithm in the two-loop two-point case, and a square root 
in the two-loop three-point case. 

By the very definition of a scalar field propagator, Pi is a quadratic form in parallel- 
space variables. Parallel-space integrations over the real line are most conveniently per- 
formed with the help of the residue theorem. To this end, we wish to close the contour 
in either the upper or the lower halfplane. This is obscured by the cut of f(Pi \ij_=k ± =o), 
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which is present in both halfplanes. The problem is readily solved by using translations 
~~ ¥ h + h an d k — > k + k\. Using the (+, — , — , — ) signature of space-time, we see that 
these translations render the propagators linear in the variables l\ and k\. Accordingly, 
the cut defined by the equation 

f(h,k 1 ) = f(P i \ l±=kx ^) = -c, (13) 

with c > 0, is now in either the lower or upper complex l\ (resp. k\) plane. This allows 
us to close the contour in the opposite halfplane. The precise location of the cut becomes 
a function of the other integration variables. This results in constraint equations for 
the contributing residues, and thus in constraints on the domains of integration for the 
remaining variables. These domains turn out to be restricted to triangles in the (/o,^o) 
plane in the case of two-loop three-point functions |27|,|30 . 



At this stage, one is left with a four-fold integral representation. There are two inte- 
grations for the remaining parallel-space variables Iq and k , and the moduli s = t]_ and 
t — k\ still have to be integrated over the positive real axis. These last two integrations 



can be carried out next. These s and t integrals have the form [^J30|1 



CO 



dt — - / ds . . ■ , . =, (14) 

t + t ±lT]Jo S + S (t) ± IT] J( at + b + irj + cs )2 _ Agt 

with the cut of the square root chosen to be along the positive real axis. The coefficients 
are functions of the parallel space variables Iq and k . With due care to the analytic 
structure of the integrand, one achieves the first integration by using either Landen 
or Euler transformations |^7|,[30]. The results allow one to interpret the t integration as 



an integral representation of dilogarithms and associated functions, and thus one obtains 
a two- fold integral representation as the final result. We once more stress that the final 
two integrations only cover finite triangular domains. They are suitable for numerical 
integrations using Gaussian quadrature. 

This programme was first successfully applied to two-loop three-point functions in 
Ref. [^TJ, and then extended to the crossed topology in Ref. [BO . For the crossed topology, 



the unavoidable presence of two propagators involving both loop momenta results in a 
more complicated cut structure. There are many more cases to be considered, but the 
calculation still follows the scheme outlined above. 



4 Results 

We now present our results for diagrams T\ through /C2 in the form of Eq. (^|). The 
number of propagators is N = 4 for N = 5 for Qi and Zj, and iV = 6 for Hi, Ji, and 
JCi. Coefficients A, B, and C are listed in Tables |I] and|^. The formulas for coefficients C 
of JF 2 , Q 5 , and TL§ are somewhat lengthy, and we have introduced the following constants: 

4 , /ox , 2, 3-V5 V5-1 , 6 ,7V 1, ,3-^5 



-_C(3) + C(2) --ln^--41n^— + — + -] + - In 
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2 
3 



•-In 3 



V5 



-In 2 



V5 



2 In' 
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T 



f, = 4 In 



3.508 941259 
2 V5-1 



_ 2v ^ln^^-5 



G, 



-1.921696 013, 

x(l — x) 
0.647466172, 



[Li 2 (— x) + lnxln(l + x)] 



o 



i/4 v^T^ 

ax — - r ln(l + x) 



Ch 



x(l — x) 
2 



7i + arctan VT5 ) + CI2 ( — 7r — arctan VT5 



v/3 



Clo 



v/5-1 



7T 



H, 



-2 In' 

2 v/3 

-0.177283 264, 
V 4 VI -4x 



In 2 



dx 

(1 — x) 2 
0.065 325 931, 



[Li 2 (— x) + lnxln(l + x)] 







i/4 VI - 4x , , 

Gtx — — m(l + x) 



CI, 



1-x) 
2 

-7r + arctan v 15 



2 

"73 

V5, 3->/5 7T /2 

In = - 

2 2 V3 V3 

-0.021441581. 



CI, 



-71 — arctan v 15 



+ 



3V3 



Ch 



In 2 + - 



(15) 



The analytic results for G5 and H5 are rather messy, and we refrain from writing them 
down here. 

In Sec. 3, we discussed the elimination of the IR-divergent pieces in Eq. (||). Substi- 
tuting in the remainder of Eq. @ the expressions for SMjj, Z w , Zjj, Ao, and vq in terms 
of Mh and A [|19j as well as the results for the various one- and two-loop diagrams, we 
may set e = and so obtain the physical result 

r- 19 

5C(2) - 3ttV3 + — + z?r(2 In 2 - 5) 



T(H -> VV) = -2\v 1 + A 



+A 5 



2C(3) + C(2) [-AZmy/l + 144 In 2 + 



1173 



4(33^ + 4^ CI 



7T 



749 



- 181tiV3 + 189L 5 - 2F 2 - 12G 5 - 18H 5 

8 

+ 12Ii + 12I 2 + 8I 3 + 72I 4 + 40I 5 + 24I 6 + 8I7 
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+ IO8J1 + 12J 2 + 12J 3 + 4J 4 + 72J 5 + 24J 6 + 24J 7 + 8J 8 + 36K! + 4K 2 

+ i U{2) (-108^^2 + 174^ - 53tt) + tt (571n2 - — - 2f 2 - 12g 5 - 18h 5 ) 

+ 40i 5 + 8i 7 + 12j 2 + 12j 3 + 4j 4 + 24j 7 + 8j 8 + 4k 2 



(16) 

where A = A/(16vr 2 ) = G F M 2 H / '(16tt 2 v^). Here, L 5 « 0.923 631827 stems from the 
all- massive lemon diagram, £ 5 = L 5 /[(47r) 4 M^] + 0(e), and is given in terms of a two- 
dimensional integral in Eq. (A86) of Ref. |TjJ, where it is called K 5 . All other constants 
are listed in Tables [I] and |2| together with Eq. (fLlp. Inserting in Eq. ([If]) the numerical 
values of these constants, we find 

T(H -> VV) « -2Aw {l + A(1.400 476 - 11.352 791 z) 

-A 2 [34.408 2(43) + 21.003 1(62) i}\ . (17) 



This agrees reasonably well with Eq. (13) of Ref. | IB|| , where the two-loop coefficient 



is found to be -[34.351(26) + 20.999(20) i]. Taking the absolute squared of Eq. (|T 
we obtain the correction factor, Ky, in the relation T(H — > VV) = KyY b{H — > VV) 
between the radiatively corrected and Born values of the H —>■ VV decay width as 

K v w 1 + 2.800 952 A + 62.030 8(86) A 2 

/ M \ 2 /A/f\ 4 
« 1 + 14.629% (^j + 16.921(2)% (p^) . (18) 

For comparison, we also list the analogous correction factor, Kf = T(H — > ff)/Ys{H 



//), for the Higgs-boson decay to a pair of fermions [13 , lj]] : 



K f « 1 + 11.058% ( I ^ 7 ) 2 - 8.908% (^) 4 . (19) 

In the derivation of Eq. (|T9|) , Yukawa couplings have been set to zero. Therefore Eq. (|T9|) 
should be a useful approximation to the full electroweak two-loop result as long as Mf 



5 Discussion and conclusions 

We are now in a position to explore the phenomeno logical consequences of our result. 
In Fig. we display the H — > VV correction factor K v of Eq. ( |TED to O(GfM^) and 
0(G 2 F Mfj) as a function of M#. For comparison, also the corresponding approximations 



for the H —> ff correction factor Kf of Eq. ( |19D are shown. In both cases, the 0{GpMjj 



terms enhance the Born results. This effect is slightly more pronounced for Ky. In the 
case of Ky, the 0{G 2 F Mj\) term acts in the same direction as the O(GfM^) term. For M H 
increasing, it rapidly gains importance relative to the OiGpM^j) term; at Mh ~ 930 GeV, 
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both terms have the same numerical value, 12.6%. In the case of Kf, the 0(G F M H ) term 
reduces and, for M H ;> 1.114 TeV, overcompensates the effect of the 0(GpM H ) term. 
Obviously, the perturbation expansions of Ky and Kf in GfM h cease to usefully converge 
for Mh values in excess of about 1 TeV. Consequently, this value may be considered as an 
upper bound on the Higgs-boson pole mass M B , if the SM is to be weakly interacting ||. 
We emphasize that this observation is independent of speculations concerning the energy 
scale up to which the SM is supposed to be valid. 

The Goldstone-boson equivalence theorem [I7| allowed us to extract power corrections 



w 



in GfM h . In the high-M# limit, these are clearly dominant over corrections in GfM ] 
and GfM 2 arising from the gauge and Yukawa sectors of the SM, respectively. For 
the H — > W + W~ and H — > ZZ decays, this was verified for the 0(GfM h ) terms in 
Refs. |3T|J3^], where the respective partial widths were completely calculated to one loop in 



the SM. There is no obvious reason why the situation should be different at the electroweak 
two-loop order. QCD corrections enter the stage only at the two-loop level, in connection 
with quark loops. The dominant two- and three- loop QCD corrections to the W + W~H 
and ZZH couplings, of 0(a s G F Mf) g§ and 0(a 2 s G F M?) |3J, respectively, are well 
under control and suppressed against the G F M H power corrections for Mh large. The 
dominance of the G F M H power corrections in the case of the ffH Yukawa couplings was 
discussed in great detail in Refs. |13| , p5| . In this context, we should mention that the 
0(GfM h ) two- loop corrections to the partial width of the loop- induced H — > 77 decay 
have been computed [|36| with the aid of the Goldstone-boson equivalence theorem. 

Since a Higgs boson with mass Mh 3> 2Mz decays almost exclusively to pairs of 
intermediate bosons and top quarks, and we have gained control over the dominant two- 
loop corrections to the respective partial widths, we are now able to predict the lifetime 
of a high-mass Higgs boson with two-loop precision. For the reader's convenience, we list 
here the Born formulas of the relevant partial decay widths: 

T B (H - W+W-) = f^|(l - r w )^ (l- rw+ \r 2 w ) , 



8vr\/2 

r.(i^zz) = ^(1-^(1-^ 

r B (H -+ ti ) = ^P^r,(l - r t fl\ (20) 

where iV c = 3 and r a = 4M 2 /M H (a = W, Z, t). We compute the total Higgs-boson decay 
width as 



H 



T B (H -> W + W~) + T B (H -> ZZ)] + K f T B (H —>tt) . (21 



We display the result as a function of M B in Fig. |3|, assuming M t = 175 GeV. Comparing 
the Born, one-loop, and two-loop approximations, the significance of radiative corrections 
to the heavy-Higgs lifetime, t b = l/r#, becomes evident. 

Equations flTSp, flTSP, and fl2"0| ) refer to the on-mass-shell renormalization scheme and 
are thus independent of the unphysical renormalization scale fi, which enters via dimen- 
sional regularization. By contrast, in other renormalization schemes such as the MS 
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scheme, there remains a weak /i dependence in the loop-corrected expressions for the 
various Higgs-boson decay widths. It is generally believed that the scheme and scale 
dependences of a calculation up to a given order indicate the size of the unknown higher- 
order contributions, i.e., they provide us with an estimate of the theoretical uncertainty. 
Applying this empirical rule to the fermionic and bosonic Higgs-boson decays, one may 
conclude that the respective perturbative results are likely to become unreliable already 
for M H £ 700 GeV |37|,|38|. In fact > for m h Z 700 GeV, the fi dependences of the MS 



results are no longer reduced if the two-loop corrections are included ||38|| . 

Finally, we should mention that Monte-Carlo studies on finite lattices provide us with 
a hint on the all-order relation between the mass and the bosonic decay width of the Higgs 



boson |39],|4(J. In Ref. [p9(1 , elastic nir scattering was analyzed nonperturbatively in the 
framework of the four- dimensional 0(4)-symmetric nonlinear a model in the broken phase, 
and the a resonance was observed. The mass m a and width of this resonance could 
be extracted from the measured scattering phases. At the same time, the pion mass m n 
and wave-function-renormalization constant Z as well as the field expectation value X in 
infinite volume were determined. In this way, it is possible to compare the nonperturbative 
value of T a with the respective perturbative value obtained by inserting m n , m a , Z, and 



£ into the tree-level formula for T a , given by Eqs. (9.4) and (9.5) of Ref. [39|. For the 
central simulation point, the nonperturbative result is found to exceed the tree- level one 
by approximately 16.7% Elf. It is interesting to confront this all-order correction with 



the one- and two-loop corrections of Eq. (|r8|). If we identify m a , Z, and £ of Ref. [^] 
with our parameters Mh, Z w , and t>o, respectively, we have 2A = (Mh/v) 2 = Z{m a /Y>) 2 , 
where the lattice constant cancels on the right-hand side. We so find that the central 
simulation point of Ref. |39| corresponds to Mh ~ 727 GeV, for which the corrections 
through 0(G F M 2 H ) and 0{G 2 F M H ) in Eq. flTJ) amount to approximately 7.7% and 12.5%, 
respectively. The tree-level, one-loop, two-loop, and nonperturbative values of K v at this 
value of Mh are indicated by the crosses in Fig. ^|. In other words, the one- and two-loop 
terms of Ky are 7.7% and 4.7%, respectively, while the higher-order terms apparently 
add up to a total of 4.3%. This is a very reasonable result, which gives support to the 
notion that the perturbative expansion in GpM H of the bosonic Higgs-decay width is 
still usefully convergent for this value of Mh, in accordance with the conclusions drawn 
above on the basis of Fig. 0. However, this observation should be taken with a grain 
of salt. The analysis of Ref. [|3£|] was performed with m n finite, while we set M w = 
for the Goldstone-boson pole mass. Nevertheless, we believe that such finite-mass effects 
will not drastically change our conclusions as long as we compare the relative corrections 
to the decay widths rather than the absolute values of the latter. This assumption is 
substantiated by the detailed inspection of the full one-loop results for the H — > W + W~ 
and H — > ZZ partial decay widths [[31~| , [32" |. 
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A Appendix: One-loop three-point diagrams 

Here, we present the relevant one-loop three-point diagrams to 0(e). We write the generic 
one-loop three-point integral as 



C(p 2 , k 2 , (p + k) 2 , m\, ml, m 2 ) 
d D q 



:2c 



(2n) D (q 2 — m\ + is) [(q + p) 2 — m\ + is] [(q + p + k) 2 — m| + is] 
The one-loop three-point diagrams in Fig. [I| are then defined as 

d = c(o,o,A4,o,ii4,o), 

C 2 = C(0,0,M 2 H ,M 2 H ,0,M 2 H ), 
^C(0,0,M 2 H ,0,M 2 H ,0), 



v 2a 
v 3 



dm 2 
_d 
dm 2 
d 



We find 
d = 



l ( 47T/i 2 

(An) 2 M 2 H \M 2 H e^ 



dm 2 

d 
dm 2 

C(2) 



2 C(0,0,M 2 H ,0,M 2 H ,0), 
2 C(0,0,M 2 H ,M 2 H ,0,M 2 H ), 



C(0,0,M 2 H ,M 2 H ,0,M 2 H ) 



iir In 2 



19 9 / ln 2 2* 

--C(3) + |c(2)ln2+i7r(-C(2) + ^ 



+ 0(e 2 ) 



V 



2d 



i ( 47T/i 2 " 

(Arr) 2 M 2 H [m 2 ^ / 
i ( A-Kfi 2 " 



;C(2) + e 



9 



C(3) 



-Cl 2 

9 V3 



+ 0(e 2 



71 

^2 +e 



S« 2 )+i—]n2 
4 2 



+ 0(e 2 



7T 



+ e 



l^ 2) + 3^3 C1 = (I 



7T 



3^ 



In 3 



(22) 



(23) 



+ 0(e 2 



(24) 



As already stated in Sec. 3, T>i and V 2b are IR divergent. We need not calculate them, 
since their prefactors in Eq. @ vanish according to Eq. (||). 
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TABLES 



Table 1: Coefficients A, B, and C to be inserted in template (||) to obtain the values of 
the divergent two- loop three-point diagrams of Fig. |l]. iV = 4 for jF i; N = 5 for Qi, and 
iV = 6 for Hi. Analytical and numerical results for F 2 , f 2 , G5, gs, H 5 , and h 5 are given in 
Eq. 
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Table 2: Coefficients C to be inserted in template (||) to obtain the values of the finite 
two- loop three-point diagrams of Fig. [I] (A = i? = 0). iV = 5 for Zj and N = 6 for J7i and 
/Cj. The numerical uncertainties in the last digit (s) are given in parentheses. 
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Q, 




Figure 1: Feynman diagrams pertinent to the H — > VV transition amplitude through 
0(G 2 F Mfj). Dashed (solid) lines represent Goldstone (Higgs) bosons. Adjacent propaga- 
tors with identical four-momenta are separated by a solid circle. 
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Figure 1: Continued. 
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Figure 1: Continued. 
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Figure 2: H -»• VV correction factor K v of Eq. (0) to 0(G F M 2 H ) and 0{G 2 F Mfj) as 
a function of M# (solid lines). For comparison, we also show the H — > // correction 
factor .fTy of Eq. ([[£]) to OiGpMfj) and 0(G 2 F Mfj) (dashed lines). The crosses indicate 
the tree-level, one-loop, two-loop, and nonperturbative values of K v at M H = 727 GeV. 
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Figure 3: Total Higgs-boson decay width T H in the Born approximation (short-dashed 
line), to 0(GfMh) (long-dashed line), and to 0(G 2 F Mjj) (solid line) as a function of M H - 
We assume M t = 175 GeV. 
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